BCS-BEC crossover in nuclear matter with the relativistic Hartree-Bogoliubov theory 
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Based on the relativistic Hartree-Bogoliubov theory, the influence of the pairing interaction 
strength on the di-neutron correlations and the crossover from superfluidity of neutron Cooper pairs 
in the 1 So channel to Bose-Einstein condensation of di-neutron pairs is systematically investigated 
in the nuclear matter. The bare nucleon-nucleon interaction Bonn-B is taken in the particle-particle 
channel with an effective factor to simulate the medium effects and take into account the possible 
ambiguity of pairing force, and the effective interaction PK1 is used in the particle-hole channel. 
If the effective factor is larger than 1.10, a di-neutron BEC state appears in the low-density limit, 
and if it is smaller than 0.85, the neutron Cooper pairs are found totally in the weak coupling BCS 
region. The reference values of several characteristic quantities which characterize the BCS-BEC 
crossover are obtained respectively from the dimensionless parameter l/(feFnffl) with a the scatter- 
ing length and fcp n the neutron Fermi momentum, the zero-momentum transfer density correlation 
function D(0) and the effective chemical potential v n - 
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I. INTRODUCTION 

The BCS-BEC crossover, which is described as the 
evolution of the pairing phenomenon from the weakly 
coupled Bardeen-Cooper-Schrieffer (BCS) type to the 
strongly correlated Bose-Einstein condensation (BEC) 
state with increasing pairing interaction, is one of the im- 
portant issues in the study of the pairing correlations in 
Fermion systems. The transition from BCS to BEC was 
investigated theoretically in several domains of physics, 
such as in excitonic semiconductors superconduc- 
tors [|J and attractive Fermi gases 0. Although BCS 
and BEC are physically two quite different limits, it was 
found that the evolution between them is smooth and 
continuous 

Since the strength of inter-atomic interaction can be 
controlled via a magnetically driven Feshbach resonance, 
ultracold Fermi atomic gases provide a special laboratory 
for studying strongly correlated Fermion systems. Re- 
cently, the BCS-BEC crossover phenomenon has been ex- 
erimentally realized in ultracold quantum atomic gas 0- 



In nuclear physics, however, the BCS-BEC crossover 
phenomenon is still a challenging topic because of the 
impossibility to control the nuclear pairing force directly. 
Indirectly, one can change the strength of pairing inter- 
action in nuclear matter by altering density or temper- 
ature. As a result, density and temperature are treated 
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as important parameters in the study of nuclear pairing 
correlations. 

For the neutron-proton pairing, the BCS-BEC 
crossover has been investigated in the 3 Si — 3 D\ channel, 
in which the strong spatial correlation and the BEC of 
the deuterons may occur at low densities f§-fllj]. 

For the neutron-neutron pairing, the correlation is ex- 
pected to be significant in low-density nuclear matter as 
well. It is well known that the bare neutron- neutron 
interaction in the So channel leads to a virtual state 
around zero energy characterized by a large negative 
scattering length a « — 18.5 ± 0.4 fm [T2|, implying a 
very strong attraction between two neutrons in the spin 
singlet state. Furthermore, theoretical predictions sug- 
gest that around 1/10 of the normal nuclear density po, 
the 1 Sq pairing gap may take a considerably larger value 
than that around po [13l4l5| . In addition, the strong di- 
neutron correlations are also supported by the enhance- 
ment of two-neutron transfer cross sections in several 
heavy nuclei (l6| . In the weakly bound neutron-rich nu- 
clei, di-neutron correlations are enhanced due to the cou- 
plings with the continuum and play an important role for 
unstable nucleus and the formation of nuclear halo 

[Ta- 
in. Recently, di-neutron emission in 16 Be was observed 
for the first time, indicating the structure of di-neutron 
clusters inside neutron- rich nuclei (2fj| . 

The progress in both theoretical and experimental in- 
vestigations on di-neutron correlations in weakly bound 
nuclei has stimulated lots of interest in searching for pos- 
sible BCS-BEC crossover of neutron pairing. The in- 
fluences on di-ncutron correlations by the Skyrme-type 
density dependent force and the finite-range Gogny force 
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are discussed in Ref. [19j. In nuclear matter, the study 
of the BCS-BEC crossover phenomenon was mainly per- 
formed by using the finite-range G ogny in teraction and 
the zero-range contact interaction [2 1| - |24| . It has been 
found that the di-neutron correlations get stronger as 
density decreases, and BCS-BEC crossover could occur 
at low densities. Besides, based on the bare force given 
by a superposition of three Gaussian functions 12 if or 
the bare nucleon-nucleon interaction Bonn potential [25[ , 
the BCS-BEC crossover in the nuclear matter has also 
been investigated. It was shown that the spatial struc- 
ture of the wave function for the neutron Cooper pairs 
evolves from BCS-type to BEC-type as density decreases. 
Erom several characteristic quantities, such as the ef- 
fective chemical potential, the quasi-particle excitation 
spectrum and the density correlation function, there is 
no evidence for the BEC state of di-neutron pairs at 
any density. However, it was argued that a Bose di- 
neutron gas may be found in the low-density limit in 
quasi-two-dimensional neutron systems as long as the 
slab is thin enough (26[. Recently, the study on the BCS- 
BEC crossover phenomenon has been extended to finite 
temperature 27 1. The thermodynamic signal from the 
temperature dependence of specific heat suggests that 
BCS-BEC crossover has been found in the low-density re- 
gion. In finite nuclei, the coexistence of BCS- and BEC- 
like spatial structures of di-neutron wave functions has 
also been revealed in the halo nucleus 11 Li [28[. From 
the di-neutron wave function, a strong correlation be- 
tween the valence neutrons appears on the surface of the 
nucleus. 

For most investigations on di-neutron correlations, 
for convenience, the phenomenological effective nuclear 
forces are used in the particle-particle (pp) channel. The 
effective interactions in the rclativistic mean-field (RMF) 
theory are also used in the pp channel [li|[3(|. However, 
one has to introduce an effective factor for this kind of 
pairing interaction in order to obtain reasonable values 
for the gap parameter [3(|. In fact, details of the effec- 
tive nucleon-nucleon interaction in nuclear matter and fi- 
nite nuclei are as yet not completely clarified, since most 
of the experimental data are not very sensitive to their 
details. It is generally believed that the bare nucleon- 
nucleon interaction should be corrected by the medium 
polarization effects ( also referred to as the screening ef- 
fects) in extracting the effective interactions in pp channel 
at various densities [3l[. Therefore, it is of great interest 
to study the behavior of di-neutron correlations and the 
BCS-BEC crossover phenomenon with different strengths 
of the pairing interaction. 

As the RMF theory has achieved lots of success in 
the descriptions of both nuclear matter and finite nu- 
clei near or far from the stability line j32Tl33 | , following a 
previous investigation [25|, the influence of the strength 
of pairing interaction on di-neutron correlations in the 
1 Sq channel will be studied with the rclativistic Hartrec- 
Bogoliubov (RHB) theory in nuclear matter in this work. 
The bare nucleon-nucleon interaction, i.e., the rclativis- 



tic Bonn potential [351 ] , will be adopted in the pp channel 
with an effective factor to simulate the medium polariza- 
tion effects and take into account the possible ambiguity 
of pairing force 3l| , and the relativistic mean- field model 
is used in the particle-hole (ph) channel. The pairing 
strength marking the crossover from the strong corre- 
lated BEC state to the weakly correlated BCS state will 
be discussed. Section [H]briefly introduces the formulism 
of the RHB theory for nuclear matter. The results are 
given and discussed in Section IIIII and finally the work is 
summarized in Section ITVl 



II. THEORETICAL FRAMEWORK 

The RMF model with nonlinear a- and w-meson self- 
coupling is employed in the ph channel to describe the 
bulk properties of the nuclear matter (32T - [3~4| | . For the 
static and uniform infinite nuclear matter, the Coulomb 
field is neglected and the space-like components as well as 
the differential of the time-like components of the meson 
fields vanish. 

In the RHB theory, the meson fields are treated dy- 
namically beyond the mean-field approximation to pro- 
vide the pairing field via the anomalous Green's func- 
tions, which gives a uniform description for the mean 
field and the pairing field [29[ . In the case of infinite nu- 
clear matter, the RHB equation is reduced to the usual 
BCS equation. The BCS ground state is defined as 



iBcs)=n( 



Uk 



v k a' kt a'_ ki ) 



(1) 



k>0 



where Uk and Vk represent the BCS variational parame- 
ters, and aL. (a_ k i) are creation operators of a particle 
with momentum fe (— k) and spin t (^) on top of the 
vacuum |— ) [36j |. 

For the ^-Sq channel, the pairing gap function A(p) 
satisfies 



A(p) = - 



1 



(27T) 3 



A • v pp (k,p) 



A(fc) 
2Ek 



dk, 



(2) 



where v pp (k 1 p) is the matrix clement of nucleon-nucleon 
interaction in the momentum space. Here an effective 
pairing interaction factor A is introduced to control the 
pairing interaction strength. The quasi-particle energy 
i?fc can be written as, 



E k = v/( £fc - M )2 + A(fc)2, 



(3) 



with the single-particle energy ek and the chemical poten- 
tial /i. The corresponding normal and anomalous density 
distribution functions are given in the following forms, 



Pk 



1 



Kk 



A(fc) 
2E k ' 



(4) 



The single-particle energy Ek follows from the RMF 
theory, 



So + Vk 2 + M* 



(5) 
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where So is the vector potential, E5 the scalar potential, 
M the nucleon mass, and M* = M + Eg the effective 



mass, 



Ss = 9a<J, S = g^o + 3 P T 3 po,3, 



(6) 



with g ai g u and g p the meson- nucleon coupling constants 
of the respective meson fields a, luq and po.3i an d T3 the 
third component of the nucleon isospin. 

With the mean-field approximation, the meson fields 
are replaced by their mean values, and could be solved 
from the corresponding equations of motion with the 
given nucleon densities, 



m 2 a 



-9aPs - ff2C 2 - g 3 a 3 , 



Wp/00,3 = 9pPb,3, 



(7) 

(8) 
(9) 



where p s , pb and pb^ are respectively the scalar density, 
vector density and isospin vector density, 



the reasonable value for the pairing gap, 



Mi 2 ) = 



q 2 + ml 



, F^q 2 ) 



A 2 



(15) 



with the meson mass and the cutoff parameter 
A<£. The vertex functions of the OBE potential r\§ and 
A<p(k,p) are listed in Table Q] For the 1 S'o pairing chan- 
nel, the matrix element v pp (k,p) is related to v pp (k,p) 
by the integral, 



pp 



(k,p) 



PP 



(16) 



where 6 is the angle between the vectors fe and p. 

In the following calculations, Bonn-B potential (3f| 
will be adopted for v pp (k,p) and the effective interac- 
tion PK1 [331 °f the relativistic mean- field theory is used 
in the ph channel. 
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p k , t k 2 dk, (10) 
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Pp = -2 / Pk.tk 2 dk, (11) 
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T 3PKl k 2 dk. (12) 



The nucleon Fermi momentum fcpi (* = is defined by 
the nucleon densities pi with = k^/Sir 2 . Therefore, for 
nuclear matter with given baryonic density pb and isospin 
asymmetry £ = (p n — Pp)/ Pb, the above equations can be 
solved by a self-consistent iteration method with no-sea 
approximation. Accordingly, the chemical potential p is 
obtained via solving Eq. ([2]), Eq. (j4|) and Eq. (fTTj) by the 
self-consistent iteration. 

The relativistic Bonn potential is used in the pp chan- 
nel, which has a proper momentum behavior determined 
by the scattering data up to high energies [3j|. It is de- 
fined as the sum of one-boson-exchange (OBE) potential 
of the six bosons <f> = a, uj, tt, p, rj, S. The vector me- 
son p includes the vector coupling channel p , the ten- 
sor coupling channel p T and the vector-tensor coupling 
channel p VT , and the pseudo-scalar mesons tt and rj in- 
clude the pseudo-scalar coupling channel it pv and rj PV , 
respectively. 

The matrix element v pp (k,p) is 



v pp {k,p) 



-A^(k,p)D lP (q 2 )F 2 (q 2 ), (13) 



where e* k is the effective single-particle energy 



:* k = \/k 2 + M* 



(14) 



D(p{q 2 ) is a meson propagator with the momentum trans- 
fer q = k —p, and F c f ) (q 2 ) is a form factor in order to get 



III. RESULTS AND DISCUSSION 

In the following discussion, the pure neutron matter 
will be mainly selected to clarify the physics. To sim- 
ulate the medium polarization effects and take into ac- 
count the possible ambiguity of pairing force [3l|, the 
effective pairing interaction factor A varying from 0.8 to 
1.4 is taken in the calculations. When the pairing gaps 
at the Fermi surface, Ap n = A(fcp n ), arc plotted as a 
function of the Fermi momentum fcp n for different effec- 
tive factor A, a maximum A(fcpn) always appears around 
fcp n = 0.8 fm _1 but changes from about 1.4 MeV for 
A = 0.8 to about 6.5 MeV for A = 1.4, and is 3.1 MeV 
for A = 1.0, i.e., the relativistic Bonn-B potential. 

For di-neutron correlations in the low density region, 
the scattering length a in the 1 Sq channel, which is de- 
fined in terms of the T-matrix for the scattering in the 
free space, is an important physical quantity. Accord- 
ing to the definition, the negative value of the scattering 
length a indicates an unbound state of neutron Cooper 
pairs while the positive value represents a bound state. 
It has been proposed to define the boundaries of the 
BCS-BEC crossover using a regularized gap equation ap- 
proach 0, Hi], [3^, |4(J , which is related with the scattering 
length a. This approach is generically applied to dilute 
systems, where the interaction matrix elements v pp (k,p) 
can be approximately treated as constant vq. Using the 
zero-energy T matrix, the relation between the scattering 
length a in the So channel and the constant interaction 
vq can be obtained as 



M 



1 



1 



dk- 



1 



Airh 2 a vo (2tt) 3 J ~'"2e(fc) : 
where e(fc) is the neutron kinetic energy, 



z{k) = s/k 2 + M* 2 - M*. 



(17) 



(18) 
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TABLE I: The vertex functions r/j, and A${k, p) of the OBE potential in the relativistic Bonn-B potential with the corresponding 
parameters gf , gf , gS , gf , ff , ff and ff . The general expression could be found in Ref. [33] - See the text for details. 



mesons 


(j> n<t> A<t,(k,p) 


scalar 


o -{gff M* J + e* k e* p - k ■ p 
S -{gf) 2 M* 2 +els;-k-p 


vector 


(gZY 2 (244 " M ") 
p v {gf) 2 2 {2els; - M* 2 ) 
p T (/ p B /2Af) 2 {M* 2 +3ele; + k-p)q 2 
p vt (f} 9 ?/M)M* q 2 


pseudo-scalar 


n PV (fn/m^Y {M* 2 + e* k £p + k ■ p) q 2 
rf v {f*/m v ) 2 {M' 2 +ele; + k-p)q 2 



The pairing gap equation ([2]) becomes 



1 = - 



(2tt) 3 



(19) 



From Eq. ([XT]) and Eq. (|19p , the regularized gap equation 
is written as 



AI 



1 



Airh 2 a 



2(2tt) 3 



dk 



1 



1 



e(k) 



(20) 



Here instead of fixing the scattering length a from its 
physical value, it is treated as a variable calculate d by 
Eq. (|2U)) for different densities similarly as in Ref. [211 ] . 
The quasi-particle energy and the neutron kinetic en- 
ergy e{k) in Eq. (|20j) arc calculated by the momentum 
dependent pairing gap A(fc) and the effective mass M*, 
which are obtained from the self-consistent RHB theory 
with Bonn potential for the pairing force. This is differ- 
ent from the regularized contact interaction model [2lj ] . 
where constant pairing gap is used. However, the differ- 
ence should be small in the low-density limit and compa- 
rable results are expected. 

In Fig. [TJ the neutron-neutron scattering length a in 
the 1 Sq channel is shown as a function of the neutron 
Fermi momentum k-p n for different effective pairing in- 
teraction factors A in the pure neutron matter. When 
A ^ 1.06, only negative branch of the scattering length 
appears and the scattering length approaches zero at high 
densities, suggesting no bound state for neutron Cooper 
pairs. However, when A > 1.08, the positive branch of 
the scattering length becomes available at the dilute den- 
sity, which implies the occurrence of a possible di-neutron 
bound state. For a given effective pairing interaction fac- 
tor of A > 1.08, with the increasing density, the scattering 
length starts from positive value and then diverges to pos- 
itive infinity, and after crossing the Fcshbach resonance, 
it reemerges from negative infinity. With increasing A, 
the Feshbach resonance is shifted to higher density. 

From the regularized gap equation (|20p . rt was 
shown [2l], HO, El| that the properties of pairing corre- 
lations can be uniquely controlled by a dimensionless 
parameter l/(fcp n a) which can give the evolution from 
BCS to BEC. The dimensionless parameter l/{kp n a) <C 
— 1 corresponds to the weak coupling BCS regime, 
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FIG. 1: (Color online) The neutron-neutron scattering length 
a in the So channel as a function of the neutron Fermi mo- 
mentum &Fn for different effective pairing interaction factors 
A in the pure neutron matter. 



while l/{kp n a) ^ 1 is related to the strong correlated 
BEC regime. From weak coupling to strong attrac- 
tion, the parameter l/(fcp n a) evolves smoothly from neg- 
ative to positive. The boundaries characterizing the 
BCS-BEC crossover can be approximately determined by 
l/(fc Fn a) = ±1 [Hill El- The unitarity limit is defined 
as l/(fcp n a) = 0, which is the midpoint of the BCS-BEC 
crossover. 

In Fig. [21 a contour plot for the dimensionless param- 
eter l/{kp n a) of the neutron Cooper pairs is shown, as 
a function of the neutron Fermi momentum fcp n and 
the effective pairing interaction factor A in the pure neu- 
tron matter. The boundaries of the BCS-BEC crossover 
{l/{kp n a) = ±1) are denoted by the dashed lines, while 
the unitarity limit (l/(fcp n a) = 0) is shown by the dot- 
ted line. Distinct features of di-neutron correlations are 
revealed for different pairing strengths. 

At low density region with fcp n < 0.2 fm -1 , the neu- 
tron Cooper pairs evolve continuously with the pairing 
strength from weak coupling BCS state to strong corre- 
lated BEC state. The effective factor A corresponding 
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FIG. 2: (Color online) A contour plot for the dimensionless 
parameter l/(fcp n a) as a function of the neutron Fermi mo- 
mentum fcp n and the effective factor A of pairing interaction 
in the pure neutron matter. The boundaries characterizing 
the BCS-BEC crossover (l/(fcp n a) = ±1) are denoted by the 
two dashed lines, while the unitarity limit between BCS and 
BEC regime (l/(kp n a) = 0) is shown by the dotted line. 

to l/(fcpn«) = 1 which characterizes the BEC bound- 
ary starts around A = 1.10 in the low-density limit and 
increases with the density. The effective factor A corre- 
sponding to l/(kp n a) = which characterizes the uni- 
tarity limit starts around 1.07 in the low-density limit 
and increases with the density monotonically. The ef- 
fective factor A corresponding to l/(fcp n a) = — 1 which 
characterizes the BCS boundary starts around A = 1.03 
in the low-density limit, decreases with the density to a 
minimum A ~ 0.85 at kp n ~ 0.35 fm _1 , then increases 
sharply. 

Besides the dimensionless parameter l/(fcpnQ)i the 
density correlation function D(q) [23|, [43|, which de- 
scribes the difference between the mean field and the 
pairing field, is also a useful quantity to study the depen- 
dence of di-neutron correlations on the strength of pairing 
force. This measure can give the transition points from 
BCS state to BEC state but not the BCS-BEC crossover 
region. At zero-momentum transfer, i.e., q = 0, the den- 
sity correlation function D(q — 0) is reduced as 

1 f°° 

D (°) = — / ( K l ~ pD k2dk > ( 21 ) 

7T Pn JO 

where Kk and pk are from Eq. ((][]). The sign change of 
D(0) has been considered as a criterion of the BCS-BEC 
crossover [23|, [43|, i.e., D(0) < means a BCS-type pair- 
ing and D(0) > represents a di-ncutron BEC state. 

In Fig. [31 the density correlation function at zero- 
momentum transfer D(0) of the neutron Cooper pairs 
is shown as a function of the neutron Fermi momentum 
fcpn and the effective pairing interaction factor A for the 
pure neutron matter. The critical line D(0) = is de- 
noted by the short dotted line in comparison with the 
reference lines of l/(kp n a) = ±1,0 obtained from Fig. [5] 




FIG. 3: (Color online) A contour plot for the density correla- 
tion function at zero- momentum transfer D(0) as a function 
of the neutron Fermi momentum fcp n and the effective fac- 
tor A of pairing interaction in the pure neutron matter. The 
boundary where D(0) = is denoted by the short dotted line. 
For comparison, the two dashed lines characterizing the BCS- 
BEC crossover and the dotted line characterizing the unitarity 
limit between BCS and BEC regime in Fig. are plotted as 
well. 

In Fig. [31 it is revealed that the density correlation 
function D(0) has a similar pattern as the dimension- 
less parameter l/(kp n a). The effective factor A corre- 
sponding to D(Q) — which characterizes the transition 
from BCS state to BEC state starts around A = 1.05 
in the low-density limit, decreases with the density to a 
minimum A ~ 1.02 at kp n ~ 0.20 fm _1 , then increases 
rapidly with the density. For A = 1.0, the results show 
that no di-neutron BEC state could occur, which agrees 
with previous results [25| . By taking the density correla- 
tion function D(0) = as a measure, the effective factor 
A characterizing the transition from BCS state to BEC 
state is smaller than the case by taking the unitarity limit 
l/(^Fnfl) = as a measure. 

In the case of the RMF theory, it has been proved 
that the neutron pair wave function in momentum space 
^pair(fc), i-e., the anomalous density Kk in Eq. satis- 
fies a Schrodinger-like equation which is expressed as (25| 

2e(fc)^ pair (fc) + 1 4 Jf fc J X-v pp {k,p)p 2 dp^ pail ( P ) 

= 2^ pair (fc), (22) 

with the corresponding energy eigenvalue 2j/ n , where v n is 
the effective neutron chemical potential obtained by de- 
ducting the momentum independent part from the chem- 
ical potential /i n , 

^ n = /i n -E -M*. (23) 

In the limit of zero density, the effective chemical po- 
tential f n behaves as half binding energy of the Cooper 
pair 0. For the evolution from the weak coupling BCS 
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FIG. 4: (Color online) A contour plot for the effective neu- 
tron chemical potential v n as a function of the neutron Fermi 
momentum fcp n and the effective factor A of pairing interac- 
tion in the pure neutron matter. The boundary of v n — is 
denoted by the short dotted line. For comparison, the two 
dashed lines characterizing the BCS-BEC crossover and the 
dotted line characterizing the unitarity limit between BCS 
and BEC regime in Fig[5]are plotted as well. 



regime to the strongly correlated BEC regime, the ef- 
fective chemical potential v n is supposed to change from 
positive to negative. 

In Fig. [U a contour plot for the effective chemical po- 
tential u a of the neutron Cooper pair is shown as a func- 
tion of /cpn and A in the pure neutron matter. The short 
dotted line denotes v n — in comparison with the di- 
mensionless parameter l/(fcp n a) = ±1,0 extracted from 

Fig.rj 

The effective factor A corresponding to v n = which 
characterizes the transition from BCS state to BEC state 
starts around A = 1.09 in the low-density limit, then 
increases monotonically with the density. For A = 1.0, 
there is no evidence for the BEC state of neutron pairs 
at any density, similarly as the case by taking density 
correlation function D(0) = as a measure. By taking 
the effective chemical potential i/„ = as a measure, the 
effective factor A characterizing the transition from BCS 
state to BEC state is larger than the case by taking the 
unitarity limit l/(k-p n a) = as a measure. 

Combining the conclusions mentioned above, it can be 
summarized that a di-neutron BEC state will appear in 
the low-density limit for A > 1.10 and there is only a 
weak coupling BCS state for A < 0.85. For A = 1.10, 
the maximum pairing gap at the Fermi surface Ap n is 
4.12 MeV around the density fcFn = 0.8 fm _1 , and for 
A = 0.85, it is 1.78 MeV. 

For symmetric nuclear matter, recent studies have 
claimed that a di-neutron BEC state can be formed at 
fcpn ~ 0.2 fm _1 after considering the medium polariza- 
tion effects [111, [23| . In order to examine the above con- 
clusion, similar investigation as the pure neutron matter 
for the symmetric nuclear matter has been done. 



For symmetric nuclear matter, by taking the dimcn- 
sionless parameter l/(fcp n a) as a measure, a di-neutron 
BEC state will occur in the low-density limit for A > 1.10 
and neutron Cooper pairs are totally in the BCS state 
for A < 0.85. By taking the density correlation func- 
tion with zero- momentum transfer D(0) as a measure, 
the di-neutron BEC state will appear in the low-density 
limit for A > 1.05. The minimum effective factor A cor- 
responding to -D(O) = is around 1.01 with the neutron 
Fermi momentum fcp n ~ 0.2 fin - . By taking the effec- 
tive chemical potential v lx = as a measure, the effective 
factor A characterizing the transition from BCS state to 
BEC state is larger than the case by taking the density 
correlation function D(0) = as a measure and the di- 
neutron BEC state may occur in the low-density limit for 
A > 1.08. Combining the conclusions mentioned above, 
for symmetric nuclear matter, a di-neutron BEC state 
will appear in the low-density limit for A > 1.10 and 
there is only a weak coupling BCS state for A < 0.85, 
which are similar as the pure neutron matter. 

In Refs. HH, HI], the BCS-BEC crossover is investi- 
gated by several characteristic quantities including the 
probability for the neutron pair partners P(d n ) as well 
as the ratios ( rms /(( n , Apn/eFn, and u n /eF n , where 
Ap n is the neutron pairing gap at the Fermi surface, 
ep n the neutron Fermi kinetic energy in Eq. (|18[) with 
CFn = e(fc = fcFn), v n the effective chemical poten- 

d n 

tial, P(d n ) = J l^'pairC^)! 2 ?' 2 ^'" with the average inter- 


— 1/3 

neutron distance d D = p n and 'l'pair(f) the neutron 
Cooper pair wave function in coordinate space, and the 
mean square radius of the neutron Cooper pairs £^ ms = 
/ Impair (^) | 2 r 4 dr/ / \^ paiiT (r)\ 2 r 2 dr. As these characteris- 
tic quantities are monotonic functions of the dimension- 
less parameter 1/ffcFna) in the regularized gap equation 
approach [2l], [40, [4l[ , they can be used to describe the 
boundaries of BCS-BEC crossover. 

In Table ITT1 are respectively listed the values of P(d n ), 
Crms/rfn, Apn/eFn and u n /eF n characterizing the BCS- 
BEC crossover in the pure neutron matter for the di- 
mensionless parameter l/(fcp n a) = 0, ±1, the zero- 
momentum transfer density correlation function D(0) 
ii. and the effective chemical potential v n — with ef- 
fective pairing factors A = 1.0, 1.1, 1.2. The correspond- 
ing values for symmetric nuclear matter are given in the 
parenthesises. For comparison, the values in the regular- 
ized contact interaction model [2l], [4(| are also given. 

In Table |TT1 it is shown that the values of P(d n ), 
Crms/^n, and y n /ep n are almost independent of the pair- 
ing interaction strength and are similar for the pure neu- 
tron matter and the symmetric nuclear matter. Fur- 
thermore, the values of P(d n ), £ rms /d n , and v n /e-p n are 
consistent with the results obtained in the regularized 
contact interaction model [2l], [4(j. In Table [TTJ the 
values of Apn/eFn slightly increase with the pairing in- 
teraction strength, which might be ascribed to the fact 
that the pairing gap Ap n increases with the pairing force 
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TABLE II: Reference values of P(d n ), £ rms /d n , AF n /eF n and Vn/eFn characterizing the BCS-BEC crossover in the pure neutron 
matter (symmetric nuclear matter) respectively for the dimensionless parameter l/(&Fn<i) = 0, ±1, the zero-momentum transfer 
density correlation function D(0) = 0, and the effective chemical potential v n — with effecti ve p airing force factor A = 
1.0, 1.1, 1.2. For comparison, the results obtained in the regularized contact interaction model [2ll. 1401] are also shown. 





A 




l/(k Fn a) 




D(0) 


V-n. 






-1 





+1 










1.0 


0.82 (0.81) 










P(d n ) 


1.1 


0.82 (0.81) 


0.99 (1.00) 




0.95 (0.95) 




1.2 


0.82 (0.81) 


0.99 (0.99) 


1.00 (1.00) 


0.96 (0.96) 


1.00(1.00) 




[21] 


0.81 


0.99 


1.00 








1.0 


0.98 (1.00) 












1.1 


0.98 (1.03) 


0.35 (0.35) 




0.52 (0.51) 




1.2 


0.98 (1.04) 


0.38 (0.38) 


0.19 (0.21) 


0.52 (0.51) 


0.26 (0.25) 




[21] 


1.10 


0.36 


0.19 








1.0 


0.25 (0.25) 










AFn/eFn 


1.1 


0.26 (0.27) 


0.74 (0.78) 




0.56 (0.55) 




1.2 


0.28 (0.26) 


0.83 (0.82) 


1.52 (1.50) 


0.59 (0.59) 


1.20 (1.21) 




[21] 


0.21 


0.69 


1.33 








1.0 


0.99 (0.98) 










^n/e Fn 


1.1 


0.98 (0.98) 


0.61 (0.61) 




0.85 (0.84) 




1.2 


0.99 (0.99) 


0.64 (0.63) 


-0.84 (-0.87) 


0.85 (0.85) 


-0.01 (-0.01) 




[40] 


0.97 


0.60 


-0.77 







strength faster than the neutron Fermi kinetic energy 
ep n does. Furthermore, difference between the present 
Afti /eF n and those in the regularized contact interaction 
model [2l| exists, which might be ascribed to the pair- 
ing force, i.e., the Bonn potential here and the constant 
interaction in the regularized contact interaction model. 
As the reference values of the characteristic quantities 
P(d n ), $ rm s/(Jn, AFn/eFn and J^n/eFn characterizing the 
boundaries of BCS-BEC crossover are obtained in the rel- 
ativistic framework by a self-consistent way and are con- 
sistent with those got by the regularized contact interac- 
tion model of non-rclativistic framework in Refs. pll. l40j . 
they provide a valuable guide in characterizing BCS-BEC 
crossover boundaries in future investigations. 

IV. CONCLUSION 

In conclusion, the influence of the pairing interaction 
strength on the di- neutron correlations in the 1 So chan- 
nel and the BCS-BEC crossover phenomenon in nuclear 
matter has been investigated based on the RHB theory. 
The effective interaction PK1 is adopted in the ph chan- 
nel and the Bonn-B potential is used in the pp channel. 
The influence of medium polarization effects on the pair- 
ing properties and the possible ambiguity of pairing force 
are simulated by an effective factor A appending on the 



Bonn-B potential. 

From the dimensionless parameter l/(k-p n ci), the zero- 
momentum transfer density correlation function D(0) 
and the effective chemical potential v ni a di- neutron BEC 
state will occur at dilute density if A > 1.10, and there 
is only a BCS state if A < 0.85. Moreover, the reference 
values of several characterized quantities P(d n ), £ rms /dn, 
AFn/eFn and f n /eF n characterizing the boundaries of 
BCS-BEC crossover are obtained in the self-consistent 
relativistic framework and are consistent with the non- 
relativistic results [2l|, Ho| , which may provide a valuable 
guide in characterizing BCS-BEC crossover boundaries 
in future investigations. 
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